Diffusion of photoexcited carriers in graphene 
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The diffusion of electron-iiole pairs, which are excited in an intrinsic graphene by the uffrashort 
focused laser pulse in mid-IR or visible spectral region, is described for the cases of peak-like or 
spread over the passive region distributions of carriers. The spatio-temporal transient optical re- 
sponse on a high-frequency probe beam appears to be strongly dependent on the regime of diffusion 
and can be used for verification of the elasic relaxation mechanism. Sign flip of the differential 
transmission coefficient takes place due to interplay of the carrier-induced contribution and weak 
dynamic conductivity of undoped graphene. 
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The electron energy spectrum of graphene has been 
characterized from marnetotransport, optical, and pho- 
toemission measurements, see reviews Q, and 
respectively. However, the relaxation and generation- 
recombination mechanisms are not verified completely at 
present time. For example, the question of a momentum 
relaxation mechanism remains controversial in spite of 
many studies, see review [1] and Refs. 5. An additional 
information on elastic relaxation can be found from opti- 
cal response under non-uniform pumping. It was demon- 
strated recently @, Q that a diffusion of photoexcited 
carriers during short time intervals can be investigated 
by the all-optical measurements using the focused pump 
and probe beams which are shifted in the spatio-temporal 
domain. A theoretical study of the charge diffusion pro- 
cess is not performed yet (only the spin diffusion, which is 
determined by a slow spin-flip scattering, is investigated 
Q) and the results of Refs. 6 and 7 remain unclear. 

In this Letter we study the spatio-temporal evolution 
of electron- hole pairs after excitation by the focused laser 
pulse. Such an evolution is governed by the energy- 
dependent diffusion of carriers caused by the elastic scat- 
tering. We analyze the case of weak intercarrier scatter- 
ing, when a peak-like distribution after mid-IR pumping 
(the energy relaxation via acoustic phonons is negligible 
up to nanosecond time scale and the carrier-carrier scat- 
tering should not be effective for low and moderate pump- 
ing levels) and a distribution spread over the passive re- 
gion, with energies less than the optical phonon energy, 
after excitation by a visible pulse takes place. The diffu- 
sion processes appear to be different for these cases and it 
is possible to extract characteristics of elastic scattering 
in low- and high-energy regions using this approach. The 
description of the optical response is performed within 
the local approach if the spatio-temporal scales of carrier 
evolution exceed the wavelength and period of radiation. 
It is important, that evolution of the transmission co- 
efficient in the high-frequency region is determined by 
the spatio-temporal evolution of the carrier energy den- 



sity because of the gapless energy spectrum with linear 
dispersion laws, see Refs. 2b and 9. Due to this, we 
have analyzed the diffusion processes for distributions of 
concentration and energy. In addition, the differential 
transmission (reflection) coefficient depends both on the 
carrier-induced contribution under investigation and on 
the weak dynamic conductivity of undoped graphene due 
to the virtual interband transitions. Such interplay leads 
to a sign flip of the relative transmission which depends 
on probe frequency and on spatio-temporal evolution of 
photoexcited carriers. Similar phenomenon was analyzed 
for the case of homogeneous geometry in the response on 
mid-IR probe. Recently, a sign flip of the spatio- 
temporal response was observed in the near-IR spectral 
region J?! but a mechanism of this peculiarity was not 
clarified. 

The spatio-temporal evolution of carriers after ultra- 
fast non-uniform photoexcitation is desctibed by the 
same distribution functions for electrons and holes /xpt ~ 
/xpt -f A/xpt. Here we took into account an effective mo- 
mentum relaxation and separated the weak asymmetric 
part of distribution, A/xpt ~ ^ (vp-Vx)/xpt/t'p where Vp 
is the momentum relaxation rate and Vp ~ vp/p is the 
velocity of carriers written through v — 10^ cm/s. The 
symmetric distribution f^pt is governed by the kinetic 
equation 



xpt 



dt 



+ (vp • Vx)A/xpt = MUt\p), 



(1) 



where the overline means the averaging over p-plane an- 
gle and the collision integrals Jr(/xt |p) describe the en- 
ergy relaxation of carriers via inelastic scattering and the 
generation-recombination processes. For the ultrafast ex- 
citation regime, when the pulse duration Tg^ is shorter 
in comparison to the diffusion and relaxation times, Eq. 
(1) should be considered Q with the initial condition at 



t ■ 



0: 



/xp,0 — 



/ 2vp — hQ 



TT / e£: 



nn 



(2) 



* Electronic address: " fedirvas@buffalo.edul 



V 7 

Here is the interband frequency of transitions due 
to the in-plane electric field E^wt exp(— irit)-|-c.c. writ- 
ten through the spatio-temporal modulated field strength 



2 



fx 1ft with the temporal form- factor wt of duration Tex- 
The Gaussian form-factor A{z) — exp(— z^)/-^?? de- 
scribes the broadening of interband transitions deter- 
mined by the phenomenological energy 27. 

Performing the averaging over p-plane in the left-hand 
side of Eq. (1), we write the kinetic equation through 
the p-dependent diffusion coefficient Dp as follows 



d 



/xpt = ^ Jr{fxt\p), 



Dp = 



2Vr, 



(3) 



Below we fit Dp through the momentum relaxation rate 
Vp determined by the model of scattering by finite- and 
short-range disorder, Up = {vdp/K)[^(jplc/K) -f uo/vd]- 
Here ^{plc/K) is a form- factor describing the finite- 
range scattering with the correlation length l^- The char- 
acteristic velocities Vd or vq correspond to finite- or short- 
range scattering contributions, respectively. Using this 
fit of Vp and the conductivity measurements, [H, Il3 we 
plot the diffusion coefficient as it is shown in Fig. 1. 
Since the density of states and Vp vanish at p — >■ 0, one 
obtains Dp (x in the low-energy region. In the high- 
energy region, at plc/^ > 1, the dependency of Dp on vp 
becomes weaker due to an interplay between finite- and 
short-range scattering mechanisms. The characteristic 
energy vh/lc is about 83 meV and 125 meV for Ic =7.5 
nm and 5 nm. [l^ The corresponding diffusion co- 
efficients are Dc ssl45 cm^/s for parameters of Ref. 10 
or Dc ~50 cm^/s and 78 cm^/s for parameters of Refs. 
11a and llb,c respectively. As a result, diffusion coeffi- 
cient exceeds 10'^ cm^/s for energies less 50 meV and this 
estimate should increase for high-mobility samples. 
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FIG. 1: (Color online) Dimensionless diffusion coefficient Dp 
in units Dc — v^lc/2va versus energy vp in units vh/lc for 
vo/vd =0.035 (1), 0.05 (2), and 0.065 (3). Dashed curve. 
Dp oc , corresponds to the short-range scattering case. 



For the case of mid-IR pumping, if hfl is less than 
the intra- and intervalley optical phonon energies hojr 
and hiJK, the distribution (2) can be used as the initial 
condition. For the case of pumping in the visible (or 
near-IR) spectral region, the initial peak (2), formed in 
the active region vp ^ Tiiok/'^, is transformed into a set 
of peaks in the passive region after a fast cascade emission 
of optical phonons. [l3| As a result, the form-factor A(z) 
in Eq. (2) should be replaced by the factor 
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Here fcp and kx are numbers of optical phonons emitted, 
EkrkK is the final energy of peak in the passive region 
corresponding to the fcr, fcif -channel, and the coefficient 
a(fcr, kfc) is determined by the ratio of emission rates for 
F and modes. If n/cur^K ^3 the number of emission 
channels exceeds 20 and the initial distribution spreads 
over the passive region. One can replace A(z) by the con- 
stant A «0.15 - 0.2 if 7 ~20 meV exceeds the interpeak 
energies in Eq. (4). 

If a diffusion time scale is shorter than the quasielas- 
tic energy relaxation scale in the passive region, vp < 
hiOK/2, then we consider the collisionless Eq. (3) with 
the initial condition (2) or /xpt=o = i^xTexA for the mid- 
IR or visible pumping cases, respectively. Using the 
Gaussian initial distribution we write the solution as 
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where Vp — ADp/l^x is the energy-dependent diffusion 
rate. For the mid-IR photoexcitation at energy hVl/2 we 
use this solution and the definitions for concentration and 
energy density 
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As a result for the Gaussian distribution of pump- 
ing intensity with the lateral size ~ lex, when £x oc 
exp[-(a;/Ze2;)^/2] and Vx = i'Re^p[-{x/lex)'^], one ob- 
tains the concentration 



rixt = 



mt + 1 



cxp 



{x /lex^ 

mt + 1 



(7) 



and the energy density Ext = fi^nxt/2 where we use 
z/ji = Vp^pn at pii = hQ/2v. For typical samples, [ll| - [l3j 



and CO2 laser pumping with HQ ~140 meV, one obtains 
the diffusion coefficient Dp^p^ =1300 - 650 cm^/s. The 
characteristic concentration nji is about 2 x 10^^ cm~^ 
under pumping intensity ^10 MW/cm^. The spatio- 
temporal evolutions of Uxt and Ext are identical for the 
monoenergetic distribution of carriers because the energy 
relaxation is omitted here. 

For the case of visible (near-IR) photoexcitation, the 
numerical integration of (6) should be performed over the 
passive region, vp < Uuk/^. Using the approximation 
(4), we rewrite the integrals (6) as follows: 
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FIG. 2: (Color online) Normalized distribution of concen- 
tration and energy density (solid and dashed curves; E = 
nvh/lc) versus dimensionless coordinate x/lex plotted for the 
time interval (0, td/2) and different correlation lengths, Ic —5 
nm (a) and 10 nm (b). 



Here n ~ 2zyflrea: A/ttZ^, the relaxation rate i^p is written 
through the first-order Bessel function of an imaginary 
argument, Ii{z), [ll[ and the dimensionless cut-off en- 
ergy is introduced as ym — ^k^c/'^u. 

Spatio-temporal evolution of concentration and energy 
density from the initial Gaussian distribution is plot- 
ted in Fig. 2 for typical correlation lengths [l^, [ij] at 
wo/fd =0.05. Quenching of energy density from initial 
Gaussian distribution at i = becomes faster in compar- 
ison to concentration if the correlation length [and ym in 
Eq. (8)] increases, c. f. Figs. 2a and 2b. The radial 
shape of the photoexcited spot transforms from Gaus- 
sian to a Lorentz-like distribution and a visible spread- 
ing takes place at time scales ^ O.Ito with r^i varying 
between 3 - 8 ps for the pumping sizes Zg^, ~300 - 500 
nm. 

The contribution of interband transitions of nonequi- 
librium carriers (described by the distribution /xpt) to 
the response on probe radiation of frequency w is de- 
scribed by the addendum to the dynamic conductivity 
(see Refs. 2b and 9; here A — ?> -fO and = huj/2v): 



^cr^tuj = — r — / dppf^pt { — - 

T^riPuj J \Puj-p + i\ 
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Plo +p + i\ 



i^E^t. (9) 



vpt 



The right-hand-side approximation is written for the 
high-frequency region, where ReAcr = and AcTxtw c>c 
uj~^ is determined through the energy density iJxt in- 
troduced by Eq. (6). It is essential that AcTxtw and 
transient optical response is determined through energy 
(not concentration) evolution. Here we used the local ap- 
proximation because /xpt varies slowly over wavelength 
distances and times ~ 27r/w. 

Finally, using the results on spatio-temporal evolution 
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FIG. 3: (Color online) Spatio-temporal evolution of differ- 
ential transmission coefficient for the mid-IR pumping case 
[energy density evolution is determined by Eq. (7)[ at probe 
energies hcj =1 eV (a), 0.78 eV (b), and 0.6 eV (c). 



of energy density given by Eq. (7) and shown in Fig. 2 we 
calculate the differential transmission coefficient for the 
above described regimes of photoexcitation. For the case 
of normal incidence of probe radiation, the transmission 
coefficient is given by 0, [13, [Tsj 
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r + ATx 
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Below we use <Jxtuj — e^/4?H-iCTtj-|- AcTxtw written through 
AfJxtw and a weak contribution to the dynamic con- 
ductivity of undoped graphene due to the virtual inter- 
band transitions, ct^. This contribution is fitted below 
as ss (e^/?i)(em/?iw — fiLo/ei) with the characteristic 
energies ~0.08 eV and ~6.75 eV. [15| Within the 
second-order approach with respect to Imcrxt(.i; the differ- 
ential transmission coefficient takes the form 



ATx 



T 



c 



2ct. 



ImA 



(l + V^ + 7re7;ic)2 



ImAcTxt^. (11) 



Both, cx E^^ and cx i?xt, contributions may be important 
here depending on an interplay between ct^ and AtTxtcj 
for the high-frequency region described by the right-hand 
side of Eq. (9). Because of an interplay between ct^ and 
AfTxiw contributions to AT/T, a sign fiip with decreasing 
of UJ is possible, see below. 

Spatio-temporal evolution of the differential transmis- 
sion coefficient is calculated below for the two regimes 
of photoexcitation with E^t given by Eqs. (7) and (8), 
respectively. We consider the single-layer graphene with 
Ic =10 nm which is placed on Si02 with -y/e ~ 1.46. In 
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FIG. 4: (Color online) Spatio-temporal evolution of differen- 
tial transmission coefficient for the pumping in visible range 
(energy density evolution is shown in Fig. 2) at probe energies 
Tiuj =1.4 eV (a), 1 eV (b), and 0.8 eV (c). 



Figs. 3a-c we plot /STy^tui/T for probe energies 1 - 0.6 eV 
and concentration ^ 2 x 10^^ cni^^. The temporal 
scale is determined by ^1.6 ps if lex ~10 /im. At 
liigli frequencies, when ct^j < 0, the ratio ATxt^j/T > 
as it is shown in Fig. 3a. A weak response with sign flip 
ATxtw/T' takes place for an intermediate frequency re- 
gion as it is shown in Fig. 3b and the ratio ATxti^/T < 
for the low frequencies, when ct^ > 0. Peak value of 
|ATxtt^/T| decreases with frequency before and after the 



sign flip region and essential spreading of response takes 
place at t > Q.b/va- Since the photoinduced concen- 
tration increases under visible excitation (we use here 
n ~ 10^"* cm~^), I ATxti^/T| also increases up to 10"'^, see 
Figs. 4a-c plotted for the spectral interval 1.4 - 0.8 eV. 
With increasing concentration the sign flip of |ATxicj/T| 
shifts to higher energies huj ^\ eV but the shape of the 
response is similar to that shown in Fig. 3. Essential 
spreading of spot takes place at i > 0.3/r/). 

The behavior of ATxtco/T at different probe frequen- 
cies under the near-IR excitation (or the differential re- 
flection coefficient IJ]) is in qualitative agreement with 
the spatio-temporal dependencies of Ref. 7. Numeri- 
cal estimates for the spatio-temporal scales in Fig. 4 
(where lex and td are about fj,m and ps, respectively; 
see above) are in good agreement with the experimental 
data. Thus, no other extra factors should be introduced 
in order to explain the sign flip peculiarities observed in 
Ref. 7. But an exact quantitative comparison requires 
a special treatment of sample parameters (e. g., doping 
level or parameters of substrate) and additional measure- 
ments including the spectral dependencies on pump and 
probe frequencies. In addition, the approximations used 
(effect of different relaxation mechanisms on distribution 
of carriers 0, [l^ and simplified description of [i3| 
) should be fitted to experimental conditions. Thus, a 
special investigation should be performed for a complete 
study of the diffusion processes under ultrafast optical 
pumping. 

Summarizing, the recent experiment and the theory 
developed here demonstrate a way for verification of the 
momentum relaxation rate in the high-energy region, up 
to the optical phonon energy. These results are also im- 
portant for the study of heating and saturation effects 
under an inhomogeneous pumping. 0, [l5| 
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